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Abstract 

The possibility to observe the medium-modified fragmentation of hard b quarks tagged 
. . . by a leading muon in ultrarelativistic heavy ion collisions is analyzed. We have found 

■ that reasonable statistics, ~ 10^ events per 1 month of LHC run with lead beams, can 

, be expected for the realistic geometrical acceptance and kinematic cuts. The numerical 

CN I estimates on the effect of the medium-induced softening b-jet fragmentation function are 

given. 
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1 Introduction 



The experimental investigation of ultrarelativistic nuclear collisions offers a unique possibility of 
studying the properties of strongly interacting matter at a high energy density. In that regime 
the hadronic matter is expected to become deconfined and a gas of asymptotically free quarks 
and gluons is formed. This is a quark-gluon plasma (QGP), in which the colour interactions 
between the partons are screened owing to collective effects [1] . One of the important tools to 
study QGP properties in heavy ion collisions is a QCD jet production. Medium-induced energy 
loss of energetic partons, the so-called jet quenching, has been proposed to be very different in 
cold nuclear matter and in QGP, resulting in many challenging observable phenomena [2]. Re- 
cent RHIC data on suppression of inclusive high-p^ charge and neutral hadron production from 
STAR [3], PHENIX [4], PHOBOS [5] and BRAHMS [6] arc in agreement with the jet quenching 
hypothesis [7] . However direct event- by-event reconstruction of jets and their characteristics is 
not available in RHIG experiments at the moment, while the assumption that the integrated 
yield of all high-p^ particles originates only from the jet fragmentation is not obvious. 

At LHG a new regime of heavy ion physics will be reached at ^/s^ = 5.5 TeV where 
hard and semi-hard QGD multi-particle production can dominate over underlying soft events. 
The initial gluon densities in Pb— Pb reactions at LHG are expected to be significantly higher 
than at RHIC, implying a stronger partonic energy loss which can be observable in various 
new channels [8, 9, 10]. In particular, the influence of the medium- modified fragmentation 
of heavy quarks on dilcpton spectra was analyzed in [11, 12, 13]. Since the estimated event 
rates for b quark production at LHG energies are expected to be high enough, in combination 
with high-pr jet production by gluon and light quark fragmentation this can give important 
information about the medium-induced effects for both light and heavy partons in nucleus- 
nucleus interactions at the LHG. 

In previous paper [14] we analyzed the possibility to observe the medium-induced softening 
jet fragmentation function ( JFF) of light partons tagged by a leading neutral or charge hadron 
in heavy ion collisions at the LHG. In this paper the possibility to measure the medium- 
modified b-JFF by a leading muon is suggested and analyzed for LHG conditions. In Sect. 2 
we give the main definitions of JFF, calculate the cross section of 5(— >■ leading fi) production 
at LHG energies with PYTHIA generator and estimate the expected event rate for the realistic 
geometrical acceptance and kinematic cuts. Sect. 3 describes shortly a model of partonic energy 
loss in QGP used to evaluate the sensitivity of b-JFF to the jet quenching effect. Discussion 
on numerical results and summary can be found in Sect. 4. 
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2 leading fi) production at LHC 

Let us recall that the jet fragmentation function D{z) determines the probability for a final 
"jet-induced" particle to carry a fraction z of the jet transverse momentum pj^. In nuclear A A 
interactions JFF for leading particles can be defined as [14] : 

^'Pt min Pt min 

where = zpj.^ = z'px is the leading particle transverse momentum, z' is the momentum 
fraction relatively to px of the parent parton (of course, without energy loss z = z' in the 
leading order of perturbative QCD), P7^*min is the minimum threshold for energy of observable 
jets. The rate of k-type jets in mid-rapidity with transverse momentum pt in A A collisions at 
the given impact parameter b is estimated as 

-^^^^^l>) - 1 # I rdrTAr,)T^(r,) ^ , (2) 

and is determined by the absolute value of partonic energy loss as well as by the angular 
radiation spectrum. Here ri^2ib,r,ip) are the distances between the nucleus centers and the 
jet production vertex y(r cos'0,rsin'?/'); rmax{b,ip) < Ra is the maximum possible transverse 
distance r from the nuclear collision axis to the V; Ra is the radius of the nucleus A; T4(ri 2) is 
the nuclear thickness function (see Ref. [15] for detailed nuclear geometry explanations). The 
effective shift Apj^{r,ip,9Q) of the jet momentum spectrum depends on the jet angular cone 
size 60 (see Fig.l). In the leading order of perturbative QCD the jet production cross section, 
oJ(jjet(k) ^ (^dp'^dy'j ^ is calculated in our case with PYTHIA6.2 [16] . The rate of high-p^- jet-induced 
hadrons is estimated as 



d{p^Ydydz' J J dpl^dy 

where the shift Apx of the hadron momentum distribution generally is not equal to the mean 
in-medium partonic energy loss due to the steep fall-off of the p^-spectrum [17] . 

For jets initiated by hght hadrons, the leading particles are the charged or neutral hadrons. 
However for heavy quark initiated jets there is the possibility to have a leading muon produced 
by semileptonic meson decays. Thus jet tagged by high-p^ muon can be identified as a heavy 
quark jet. Note that ~ 20% of i?— mesons and ^ 12% of D— mesons decay to muons, about 
half of the muons from i?— decays being produced through an intermediate D [18]. 
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We used PYTHIA6.2 [16] with CTEQ5L pdf parameterization to calculate the cross section 
of b-jet production and the corresponding spectra at ^/s^ — 5.5 TeV and to estimate the 
expected event rate for the realistic geometrical acceptance and kinematic cuts. To be specific, 
the geometry of Compact Muon Solenoid (CMS) detector is considered [19, 20]: the pseudo- 
rapidity coverage \ rj \< 3 for jets and | |< 2.4 for muons. We define the muon as a 
leading particle if it belongs a hard jet and carries larger 20% of the jet transverse momentum. 
To be specific, the jet energy is determined here as the total transverse energy of the final 



particles collected around the direction of a leading particle inside the cone R = \/Arp~+"A^ — 
0.5, where rj and </? are the pseudorapidity and the azimuthal angle respectively. Extra cuts 
Pj. > 5 GeV/c and E'^^* > 50 GeV were apphed. Then the corresponding pp cross section for 
S(— > leading //) production is ^ 0.7 pb, and Pb— Pb cross section is estimated as 0.7 pb x 
(207)^ fa 0.03 mb. The corresponding event rate in a one month Pb— Pb run (assuming 15 days 
of data taking), if = 1.3 x 10^ s, with luminosity L = 5 x 10^^ cm~^s~^, is A^ev = HapbPbL ^ 

2 X 10^ in this case. Increasing the minimal jet energy results in reducing expected statistics, 
e.g. for Eq^^ > 100 GeV the estimated rate is only ~ 10^ events. 

3 The model for simulation of jet quenching 

In order to test the sensitivity of b-JFF to the jet quenching, the following event-by-cvent 
Monte Carlo simulation procedure was applied (see Refs. [15, 21] for details of the model). 

• Generation of the initial parton spectra with PYTHIA (fragmentation ojf). 

• Generation of the jet production vertex at the impact parameter b according to the 
distribution 

dN^ ^ T4n)TA{r2) 

rlihrlr ^tt fmax ' ^ ' 

"^"^"^^ Jdij J rdrTA{n)TA(r2) 



• Calculation of the cross section a — J dtda/dt for scattering of a parton with energy E off 
the "thermal" partons with energy (or effective mass) mo 3T <^ E {T is the medium 
temperature) and generation of the transverse momentum transfer ti according to the 
distribution 

da ^ 27rag(t) E' Utt 

dt E^-my ' {33-2Nf)ln{t/AlcD) ' ^ ^ 

Here C = 9/4, 1,4/9 for gg, gq and qq scatterings respectively, is the QCD running 
coupling constant for Nf active quark flavors, and Aqcd is the QCD scale parameter which 
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is of the order of the critical temperature, A-qcd Tc ~ 200 MeV. The integrated cross 
section a is regularized by the Debye screening mass squared IJ-d{T) ~ 47rQ;sT^(l + A^//6). 
The maximum momentum transfer imax = ~ ("^9 + "^■o)^][s ~ ("^g ~ nT'oYj/s where 
s — 2moE + ml + rri^, rUq is the hard parton mass. 

• Generation of the transverse distance between scatterings, li — (rj+i — Tijpx/E: 

— = A-^(t,+i) exp (- I X-\n + s)ds) , X-\r) = a{r)p{r) , (6) 

where r is the proper time, A is the in-medium mean free path, p oc is the medium 
density. 

• Reducing the parton energy by coUisional and radiative loss per scattering i: 

A-Etot,i = A£^col,i + AE^rad.i , (7) 

where the coUisional part is calculated in the high-momentum transfer approximation, 

and the radiative part is generated according to the energy spectrum dl/dcu obtained in 
the frame of BDMS model [22] generalized to the case of heavy quarks - the "dead cone" 
approximation [23] (but note there are exist more recent developments on heavy quark 
energy loss in the literature [24, 25]): 

— \ 1 dl 1 - llX'^ (^V'^ 

dw " (1 + {lujf/^Y du '"''=° ' \^xl) \e) ' 



duj ttLo; 2 



In [cos (cuiTi)l , (10) 



where ti = L/{2Xg), y — uj/E is the fraction of the hard parton energy carried by the 
radiated gluon, and Cr — 4/3 is the quark color factor. A similar expression for the 
gluon jet can be obtained by substituting Cr — 3 and a proper change of the factor in 
the square bracket in (10), see Ref. [22]. The allowed range of values Ui — AE'rad.i in (10) 
is from cUmin = Elpm = jJ^'oXg, the minimal radiated gluon energy in the coherent LPM 
regime, to initial jet energy E. 

• Calculation of the parton transverse momentum kick due to elastic scattering i: 

A^. = {E- -^f - (p - - ^1)2 _ ml (12) 
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• Formation of the additional (in-medium emitted) gluon with the energy uji — AE'rad,? and 
the direction relatively to the parent parton determined according to one of two possible 
simple parameterizations for the emission angle 9: the "small- angular" parameterization, 



where 9o ~ 5° is the typical angle of the coherent gluon radiation estimated in [26]; or 
the "wide-angular" parameterization, 



• Halting the parton rescattering if a parton escapes from the dense zone, or 2) QGP 
cools down to = 200 MeV, or 3) a parton loses so much energy that its PTi^) drops 
below 2T{r). 

• In the end of each event adding new (in-medium emitted) gluons into PYTHIA parton 
list and rearrangements of partons to update string formation are performed. 

• Formation of the final particles by PYTHIA (fragmentation on). 

The medium was treated as a boost-invariant longitudinally expanding quark-gluon fluid, 
and partons as being produced on a hyper-surface of equal proper times r [27]. In order to 
simplify numerical calculations in the original version of the model we omit the transverse 
expansion and viscosity of the fluid using the well-known scaling solution due to Bjorken [27] 
for a temperature and density of QGP at T > Tc ~ 200 MeV: 



For certainty we used the initial conditions for the gluon-dominated plasma formation expected 
for central Pb-Pb coUisions at LHC [28]: 

To ~ 0.1 fm/c. To ~ 1 GeV, Pg ^ 1.957^ . 




(13) 



airy/' = SoTq/', T(r)r^/3 = Tot'^^', p{t)t = p^r^. 



(15) 
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4 Numerical results and conclusions 



Figure 2 shows the b-JFF (1) of leading muons for the cases without and with medium-induced 
energy loss in central Pb— Pb collisions with the two parameterizations of the distribution on 
the gluon emission angles (13) and (14); the same geometrical acceptance and kinematic cuts 
as described in Sect. 2 were used. One can see the softening b-JFF due to partonic energy loss 
at 2; > 0.4. The effect enhances with z decreasing (see Fig. 3) and is more pronounced for the 
small-angular radiation. The reason for the latter fact is following. The contribution of the 
small-angular radiation to the total jet energy loss (due to "out-of-cone" partonic energy loss) 
is much lesser as compared with the broad-angular radiation. The former docs not disappear 
totally mostly because not only leading (parent) parton, but all partons of a jet pass through 
the dense medium and emit gluons under the angles Q relatively to their proper directions, 
which in general may not coincide with the jet axis (determined by the direction of a leading 
particle) and sometimes be even at the jet periphery. The broad-angular radiation increases the 
"out-of-cone" part of partonic energy loss and thus decreases the final jet transverse momentum 
pJ^* (which is the denominator in the definition oi z= Pt/Pt J^F (1)) without any influence 
on the numerator of z and, as a consequence, in reducing effect on JFF softening. 

Note that in the real experimental situation the jet observables will be sensitive to the 
accuracy of jet energy reconstruction in a high multiplicity environment, in particular, to the 
systematic jet energy loss. However, since the average reconstructed jet energy in Pb— Pb 
collisions is expected to be the same as in pp interactions (see section "Jet detection at CMS" 
in Ref. [8]), the short measure of jet energy will be the well-controlled systematic error for heavy 
ion as well for pp collisions, and it can be taken into account using the standard calibration 
procedure. 

In summary, the channel with the muon tagged b-jet production in ultrarelativistic heavy 
ion collisions was first analyzed. The reasonable statistics, ~ 10^ events per 1 month of LHC run 
with lead beams, can be expected for the realistic geometrical acceptance and kinematic cuts. 
The effect on the medium- modified b-jet fragmentation was numerically studied for Pb— Pb 
collisions at the LHC. The significant softening b-jet fragmentation function determined by the 
absolute value of partonic energy loss and the angular radiation spectrum is predicted. 

Acknowledgments. Discussions with Yu.L. Dokshitzer, O.L. Kodolova, C. Roland, I.N. Var- 
danyan, R. Vogt, B. Wyslouch, B.C. Zakharov and CM. Zinovjev are gratefully acknowledged. 
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Figure 1: The schematic view of a jet with gluons emitted inside and outside jet cone. 
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Figure 2: B-jet fragmentation function for leading muons without (solid curve) and with medium- 
induced partonic energy loss for the "small-angular" (13) (dotted curve) and the "broad-angular" (14) 
(dashed curve) parameterizations of emitted gluon spectrum in central Pb— Pb collisions. Applied 
kinematic cuts are described in the text. 
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Figure 3: The ratio of B-jct fragmentation function for leading muons with energy loss to one without 
energy loss in central Pb— Pb collisions. The dotted curve is the result for the "small-angular" radia- 
tion, the dashed curve - for the "broad-angular" radiation. Applied kinematics cuts are described in 
the text. 
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